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The Roothaan and Pople-Nesbet approaches for real atoms are adapted to quantum dots in 
the presence of a magnetic field. Single-particle Gaussian basis sets are constructed, for each dot 
radius, under the condition of maximum overlap with the exact functions. The chemical potential, 
the charging energy and the total spin expected values have been calculated, and we have verified 
the validity of the quantum dot energy shell structure as well as the Hund rule for electronic 
occupation at zero magnetic field. For finite field, we have observed the violation of Hund's rule 
and studied the influence of magnetic field on the closed and open energy shell configurations. 
We have also compared the present results with those obtained with LS-coupling scheme for low 
electronic occupation numbers. We focus only on ground state properties and consider quantum 
dots populated up to 40 electrons, constructed by GaAs or InSb semiconductors. 

I. INTRODUCTION 

The influence of the spatial confinement on the physical properties, such as electronic spectra of nanostructured 
systems, is a topic of growing interest. Among the several kinds of confined systems one can detach low dimensional 
electronic gases and impurity atoms in metallic or semiconductor mesoscopic structures^ as also atoms, ions and 
molecules trapped in microscopic cavities ^Sii 4 ^ where the effects of the confinement become important when the 
typical quantum sizes - such as Fermi wavelength - reach the same order of magnitude as the sizes of the cavities. 
However, the energy spectrum of such systems is not only determined by the spatial confinement and the geometric 
shape, but also by environmental factors such as electric and magnetic fields. It is also defined by many-body effects 
as electron-electron interaction, which may even become more important than the confinement itself. In all cases, a 
correct description of the physical properties of the problem requires the system wavefunction to reflect, in an appro- 
priated way, the presence of both confinement, internal and external interactions, and the corresponding boundary 
conditions. Wc should also mention the existence of other confined systems, for example, the ones constituted by 
phononspi plasmonsi or confined bosonic gases£ 

Low dimensional electronic gases are defined in semiconductor structures when the bulk translation symmetry is 
broken in one or more spatial dimensions, giving origin to 2D (quantum wells), ID (quantum wires) or 0D (quantum 
dots) systems. In such structures, the charged carriers loose, for some interval of energy, the characteristic of being 
delocalized in both three spatial dimensions and become confined in regions of mesoscopic sizes inside the crystal. 
This fact transforms the continua energy bands into broken subbands and energy-gaps, or even into fully discrete 
energy states as for example in semiconductor quantum dots (QDs), which are the main type of confined system to 
be addressed in this work. 

An important point to have in mind is that the usual charging modelfSiiSiiiiiSi which reduces the electron-electron 
interaction inside the dot to a constant proportional to its electronic occupation, is able to reasonably reproduce the 
experimental findings for metallic dots. However, in order to obtain a more realistic description of many-particle semi- 
conductor dots, their relatively lower electronic density makes imperative to consider the electron-electron interaction 
under a microscopic point of view. In strong confinement regimes, it may even be included in a perturbative scheme. 

Various are the approaches that have been used to deal with many-particle QDs. Among them, one can cite charging 
model, correlated electron model^ Green functions^ Lanczos algorithm^ Monte Carlo method?^ Hartree-Fock 
calculationsi2ii2*i2i2£ and density functional theoryiSi It is useful to emphasize that not only the many-body effects 
but also the spatial symmetry (geometric confinement) become indispensable ingredients for a precise determination 
of quantum effects in the electronic structure of semiconductor QDs. 

Among the several kinds of geometries confining a QD, maybe the most common be a two-dimensional one defined 
by a parabolic potentialiSSiSiS 4 * 2 ^ 2 ^ Here we will consider a three-dimensional QD defined by an infinite spherical 
potential. The former describes QDs lithographically defined in the plane of a two-dimensional electron gas, while the 
latter describes QDs grown inside glass matrices. Some of the commonly studied topics in 3D QDs are the formation 
of energy shells in their spectra^ the control of electronic correlations^ the formation of Wigner molecules in the 
system jS, and the influence of the Coulomb interaction in their spectra^S^i In such spherically defined 3D QD, both 
spin and orbital angular momenta are good quantum numbers, and the many-particle eigenstates can be labelled 
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according to the usual LS-coupling scheme^ where the exact analytic many-particle eigenstates are given as a sum 
over appropriated Slater determinants. However, in order to deal with QDs having higher occupation, the LS- 
coupling is no longer appropriated. Therefore, we have chosen to use in this paper the Roothaan and Pople-Nesbet 
matrix formulations 33 of the single determinant self-consistent Hartree-Fock formalism, respectively appropriated for 
handling closed and open shell configurations. With them, we also calculate the QD addition spectrum and show how 
a magnetic field is able to violate the Hund rule. Here we consider a QD populated with up to N — 40 electrons under 
the presence of a magnetic field, and expand our basis in a set of properly optimized Gaussian functions. As a note 
one could, in principle, also include spin-orbit coupling in the model, where L and S would then no longer be good 
quantum numbers, but for an infinite spherical potential such coupling yields no contribution to the total energy of 
the system. 

The paper is organized as follows. In section II, for completeness and in order to introduce our notation, both 
restricted and unrestricted formalisms are resumed. In section III we show the Hamiltonian model and discuss how 
both QD chemical potential and charging energy are calculated; also, we present details on how the inclusion of a 
magnetic field changes the previous formalisms, as well as how the Gaussian basis set is constructed. Finally, in 
section IV we show our results and we use section V to line up our conclusions. 



II. THEORETICAL METHOD 



The Hartree-Fock (HF) approach assumes that the TV-electron ground state of an interacting system is given by the 
single Slater determinant \^q) — \XiX2---XaXb---XN— iXn)i where the set of optimized spin orbitals {x a l a = l-..iV} is 
obtained by the minimization of the total energy Eq of such state, 
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where the kinetic h aa , direct J a b, and exchange K a b contributions to Eq are given by 

h aa = (Xa\ K \ Xa ) = ~ J d^ lX * a (l) [V 2 Xa (l)] , 



Jab = ( Xa \ Jb \Xa) = ~ d<lx lXa(l) 
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with m (e) being the effective electron mass (material dielectric constant), and V 2 = r~ 2 (<9 / dr(r 2 d / dr) — L 2 ) being 
the Laplacian. The symbol d 4 x is used in the integrals because the orbitals {x a } involve both spatial and spin parts. 
The minimization of Eq. (JTJ yields the self-consistent integro-differential HF equation, whose procedure in order to 
obtain the optimized spin orbitals can be found elsewhere^ it is given by / \x a ) — £ a lx a ); where the Fock operator 
is 



/ = K + Y, t J b - K b\ 



(3) 



It can be straightly verified that the application of / on \\ a ) can be obtained from Eq. J2J). The ground state energy 
Eq is calculated in each iteration of our code until a convergence of 10~ 9 eV be reached. 

The above expressions are general in the sense that no supposition was taken on the spin orbitals {x a }- Furthermore, 
the confinement potential V(r) should be added to H, however, since we have chosen to work with an infinite spherical 
potential barrier, we have excluded it from our expressions. 



A. Roothaan approach 

The first supposition, which is the core of the Restricted Hartree-Fock approach (RHF) , consists in assuming that 
both a (up) and (3 (down) spin functions are restricted to have the same spatial function or, more specifically, Xi( x ) = 

{ipj(r)a(Lj), ifjj(r)P(uj)} . In such way, the ground state of the system becomes \^q HF ) — tp 1 Tp 1 ...^ a Tp a ...^ N / 2 7 pN/^ 

where a spatial function with (without) upper bar labels the spin-down (spin-up) state. Therefore, within the RHF 
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approach, the spatial orbitals {ip a \a = 1...N/2} are doubly occupied and HF ) represents a closed shell system 
having total angular momenta J = L = S = 0. After integration over the spin degree of freedom^ the HF equation 
in the restricted formalism becomes / \ip a ) — £ a \ip a )i where the Fock operator is given by 

N/2 

f = K + I 2 ^ - Kb] , (4) 

b 



and the ground state energy becomes 



N/2 N/2 N/2 

E* HF =2Y d h aa + Y d Y j \2J ah -K ab ]. (5) 

a a b 

Here, h aa = {tp a \h a \ip a ), J ab = (ip a \ J b \ip a ), and K ab = (ip a \ K b \ip a ); so, the integrations in Eq. J5J involve only 
spatial functions in the RHF approach. 

The Roothaan formalism consists in translating the RHF equations into a matrix formulation. This is achieved by 
expanding the set {tpi\i = l.-.fc} to be minimized, in a set of known basis functions {4> v \v — l...fc}, 

ipi = ^2C vi (f> I/ , (6) 

V 

where the spatial orbitals having k > N/2 are empty. Within this procedure, the coefficients C v % become the 
parameters to be iterated. When Eq. (JSJ is inserted in the closed shell HF equation (/ given in Eq. Q) and 
multiplied from the left by (p* (r), one obtains the Roothaan characteristic huh matrix expression, 

FC = SCe, (7) 

where S is the positive defined overlap matrix between the basis functions, whose elements are 

V = J d 3 r^(r)^(r), (8) 

C is the matrix of the expansion coefficients CV;, whose columns describe each spatial orbital ip^, e is the diagonal 
matrix of the orbital energies £,-, and F is the matrix of the Fock operator /, whose elements are 

= J d 3 r^(r)/^(r). (9) 

In order to explicitly write the Fock matrix, it is convenient to introduce the charge density of the system. For a 
closed shell configuration, described by a single determinant where each spatial orbital is doubly occupied, it can be 
written as p(r) — 2^^ 2 |?/> a (r)| 2 = J2u ^/^^v( r )^C( r )! where one defines the density matrix P whose elements 
have the form 

N/2 

Pfiv — 2^C Aia C* a ; (10) 

a 

observe that an integration of p(r) over all space yields N. The elements of the Fock matrix are obtained when one 
inserts Eq. J3J into Eq. I© and uses Eqs. © and ljTU|) . After some algebra, one gets = T^ v + G MI/ , where the 
kinetic contribution is given by 



(11) 



while the Coulomb contribution is written as 

2 



A a 



(12) 



-\ J d 3 n J d3 r2 ^ (ri) ^ (r2) ^_^^ (ri) ^ (l . 2) 
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The self-consistency of this method appears as the dependence of F on P that, in turns, depends on C, which are the 
parameters to be determined. 

Therefore, the procedure for the solution of Eq. (JJJ) must follow the steps: i) Given a confinement potential for the 
system, one specifies N and {<f>^}', h) The integrations in and are done; iii) An initial guess is used for P; iv) 
With P and the two-electron integrals one obtains the matrix G, which is added to the matrix T (once calculated, 
this one-electron matrix no longer changes, since it does not depend on P or C) to form F; v) F is diagonalized in 
order to get C and £, and the solution is used to form a new matrix P from Eq. I|1(J|I : vi) This iteration is performed 
until the desired convergence be found. The criterion for convergence may also be set in the change of the ground 
state energy that is calculated in each iteration; in the Roothaan formalism it can be written as 

E o HF ^EE^ + F »») ■ ( 13 ) 



B. Pople-Nesbet approach 

The second supposition is based on the relaxation of Roothaan's restriction by letting the a and (3 spin functions 
have different spatial components, which is the core of the Unrestricted Hartree-Fock approach (UHF). Thus, Xi( x ) = 

yip" (r)a(w), ipj (r)/3(w)|, where spin-up (spin-down) electrons are described by the spatial orbitals (Cb = 

{{^jlj — l—k}). Within the UHF formalism, an unrestricted wavefunction has the form ^ UHF ^ = tp"ip^...^, 

which represents an open shell system since no spatial orbital can be doubly occupied. These UHF functions are not 
necessarily eigenstates of the system having defined L and 5* values, and N is no longer restricted to an even number, 
but it has to satisfy the condition N = N a + N@ , the sum of spin- up and spin-down electrons. 

In analogy with the RHF case, the spin can be integrated in the UHF formalism^ The main difference here is 

that there are two coupled HF equations to be simultaneously solved. They are given by f a ^ 3 ip^j — /) 

where the respective Fock operators are 

r" = ^ + £ [j a J p + £ J a /a - (14) 

a a 

Notice that this equation reduces to Eq. if ijf? — <^ The way in which the corresponding terms in f a /P operate 

on ^"^^ can be obtained from Eq. (J2J, keeping in mind that here one deals with only spatial contributions. In 

both f a and f 13 there is presence of a kinetic term hj, a direct J« and an exchange K% term between electrons 
with same spin, and a direct term Ja^° between electrons with different spin. Observe that it is such interdependence 
among f a and ipj, as well as among f@ and ip", that makes necessary a simultaneous solution of the two HF equations 

in the UHF approach. Such solution yields the sets {ipj } and |^ | that minimize the ground state energy Eq HF of 
the unrestricted state \^q HF "), given by 

e u » hf = £ K a + £ *L + ~ £ £ w - K"} + \ E E VI - «%] + E E a 5 ) 

a a a b a b a b 

Here, h a J/ = (C //3 | K |c //3 ), J a J = (CI 4 IC) = (<| ^ \<) = C -4T = (CI J° IC> = (CI J? IC> = 
J£f, plus the analog term for J%. Also, Kg = (CI K% |C> = (C I IC) = Kg, P^s the analog term for K%; 
observe that there is no K"^ term since such interaction does not exist. Furthermore, the explicit expression for each 
one of these integrals is straightly obtained from the spatial contributions given in Eq. J5J. As should be expected, 
Eq. O reduces to Eq. © when ip" = if)?. 

The Pople-Nesbet formalism has analogy to the Roothaan one. The difference is that, instead of Eq. 1(5} , there are 
two distinct expansions for the spatial orbitals C and C in terms of the known basis functions {(p^v — l...fc}, given 
by 

C //3 = £C/V as) 
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By following the same procedure as in the Roothaan matrix expressions, one obtains the two characteristic k x k 
coupled equations for the Pople-Nesbet approach, 



pa//3(-<a//3 _ gQa / '/3 £ a / '0 _ 



(17) 



Here the overlap matrix S is defined in Eq. (jHJ, and the meanings of C a ^ , e a ^, and F a ^ become understood from 
the discussion following that equation. The respective Fock matrices are given by 

Kl P = / d 3 <(r)r/^(r). 



(18) 

It is also convenient in this UHF formalism for open shell systems to introduce the charge density of the up and down 

spin electrons, p a ^(r) — J^a ' i 3 a^i r ) ~ EuEi/^'^W^W' where the elements of the corresponding 
density matrix for the up and down states are 



jV' 



</f> 



pa/0 _ na/fina/ 



10* 



(19) 



From these expressions one can define two new quantities. One is the total charge density, p T (r) — p a (r)+p^(r), which 
yields the total number of carriers, N, after integrated over all space; the other is the spin density, p s (r) = p a (r)— p^(r), 
whose integral over all space yields 2Ms- The latter shows that unrestricted wavefunctions are eigenfunctions of Sz, 
but not necessarily of S 2 . Consequently, one can define the total, P T , and spin, P s , density matrices for the system 
as 

pT = pa + pP_ 

P s ^P a P> 3 . (20) 

Furthermore, the elements of the two Fock matrices can be obtained by inserting Eq. (|14H into Eq. I|18|) , and using 
Eqs. H16|) . (|19|) and J20J). After some algebra, one finds Fjfv = T^ v + G^if , where the kinetic matrix is defined 
in Eq. Ijlll) . and the two Coulomb matrices G"(f are given by 

G % P = 7 EE \ P Lj d ^ J rf 3 r 2 ^(r 1 )C(r 2 )^-i^^(r 1 )0 A (r 2 ) 



-Pxi J rf3 n / d 3 r 2 ^(n)C(r 2 )— l_^ A ( ri )^(r 2 ) 



(21) 



Notice that the self-consistency again lies in the fact that both F matrices depend on the P matrices which, in turns, 
depend on C matrices. Finally, the coupling between up and down spin equations appears on the property that F a 
(F@) depends on V fi (P Q ) through P T . 

The procedure for the solution of Eq. (|17|) is the same used for Eq. J7J, having in mind that now a simultaneous 
solution must be found for the two distinct spatial orbital sets. If the criterion for convergence in the Pople-Nesbet 
formalism is the ground state energy, it can be written as 



pUHF _ I \pT rp , pa pa , p/3 p(3 ] 



(22) 



As already mentioned, one disadvantage of this approach is that unrestricted functions, in general, are not eigen- 
statcs of the total spin S, even though they have defined Sz values. However, one may use the expressions^ 3 . 



<* 2 > 



UHF 



N a -N l3 \ f N a -N (i 



N a N 13 



1 +^-EE 



/ < / 4 fia ^ v b°^ 



in order to get an estimative for S, while Sz is obtained from 

(Sz)uHF = 9 E E _ ^») 



(23) 



(24) 



As a last observation, the expected value of S will, in general, be higher because of the contamination from other 
symmetries. 
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TABLE I: Values of the five coefficients and exponents that optimize the set of Gaussians used in the expansions for the six 
orbitals considered in this work. 



orbital 


Vi, D 1 (1Q- 4 ) 


V 2 , D 2 (10- 4 ) 


V 3 , D 3 (10~ 4 ) 


Vi, £> 4 (10- 4 ) 


Vi, £> 5 (10- 4 ) 


Is 


0.3229, 0.003 


0.2353, 0.328 


0.4554, 0.043 


0.9069, 0.330 


0.8317, 0.006 


lp 


0.5902, 0.762 


0.0646, 0.564 


0.5941, 0.447 


0.4599, 0.350 


0.9727, 0.724 


Id 


0.7301, 1.120 


0.9570, 0.915 


0.8227, 1.620 


0.1950, 0.492 


0.4194, 1.100 


2s 


0.4823, -2.260 


0.9875, -2.440 


0.1082, -1.530 


0.4006, -2.840 


0.0342, -2.260 


If 


0.6377, 1.500 


0.4539, 1.920 


0.9730, 1.440 


0.4921, 1.610 


0.5498, 1.360 


2p 


0.8513, -1.420 


0.0413, -3.720 


0.0565, -0.401 


0.8797, -1.410 


0.0457, -1.700 



III. SPHERICAL QUANTUM DOT 



As an application of both RHF and UHF formalisms, we consider a QD with radius Rq confined by an infinite 
spherical potential in the presence of a magnetic field B = £?o(0,0, 1) that can be populated up to 40 electrons. The 
single-particle Hamiltonian has the form 



H « = f ( k + f A ) 2 + -?ir B - s > 

2m \ nc / h 



(25) 



where fi B = eh/(2moc) is the Bohr magneton, g is the bulk g- factor, and the vector potential is used in the symmetric 
gauge A = (B x r)/2. By using the atomic units, Er v = e 2 /(2ao) for energy and a — h 2 /(m e 2 ) for length, the 
Hamiltonian Hq can be written as 



„ _ 1 o§ 
m Rq 



)_d_ ( 2 d_\ L2 
' x 2 dx \ X dx) + x 2 + 21% { 2 



mgSz) 



|Msin 2 (0) 



(26) 



where m = m/nio, Ib — y hc/(eBo) is the magnetic length, and x — r/Ro is a dimensionless variable. Without 
magnetic field, the normalized spatial eigenfunctions of Hq are given by 



U r ) = <PnJ,mi( x > 



1 



1 1/2 



[ji+i (a n i)Y 



ji (otnix) Yl >mi (9, ( 



(27) 



The boundary condition at the surface r = Ro determines a n i as the n th zero of the spherical Bessel function ji(a n ix); 
also, the spherical harmonic Yi m . (6, <f>) is the well known eigenstate of L 2 and Lz- 



The Hamiltonian for the electron-electron interaction, in units of Er v , has the form 



_ 2 a 

■Ubk — —- 



1 



e Ro Ixi - x 2 



(28) 



By using the multipole expansion, 



4tt 



Xi - x 2 



^ 2k 

K = 



E (- 



(29) 



all the angular part of the problem can be solved analytically and, after inserted into our numerical code, we are left 
to solve only the radial degree of freedom. 

Without magnetic field, we take into account in Eqs. © and l |16(l the spatial orbitals that define the six lowest 
energy shells (Is, lp, Id, 2s, If, 2p) under this symmetry^ Therefore, the index v = n,l, mi can assume up to 40 (20 
spin-up and 20 spin-down) possible values for the states within those shells. Certainly, the inclusion of a magnetic field 
lifts both spin and orbital degeneracies of those states. In our application, we shall consider two possible materials 
forming the QD, namely GaAs (wide-gap material) and InSb (narrow-gap material), whose defining parameters are 
m = 0.065, g = 0.45 and e — 12.65 for GaAs, and fh = 0.013, g = —53.1 and e — 16.5 for InSb. 

Due to the presence of a magnetic field, some modifications on both RHF and UHF expressions must be made. 
In the Roothaan approach, the quadratic term in Bq (~ l B 4 ) must be added to the single-particle contribution 
present in h aa in Eq. J5J, and to the matrix T^ v in Eq. JUJ; the linear terms in Bo (~ l B 2 ) are zero for a closed shell 
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FIG. 1: Comparison between the expanded Gaussian basis sets used in our calculations and the respective exact spherical 
Bessel functions, for a GaAs QD having N = 2 and Ro = 100 A. The left panel shows the probability densities, while the right 
one shows the wavefunctions. The agreement is good for all orbitals except very closed to the origin for s-states. 

configuration, thus yielding no contribution. However, all terms proportional to Bo have to be considered in the Pople- 
Nesbet formalism, since we will be dealing with an open shell configuration. Both linear (~ BqLz) and quadratic 
(diamagnetic) terms can be straightly added to the definitions of h aa (Eq. (J2J) and (Eq. Moreover, the 

inclusion of the spin-dependent linear term (~ BqSz) in h aa will impose the kinetic matrix be decomposed in its 

components T^l^ ', as it occurs with G^J 3 in Eq. (|21|l . Thus, under a magnetic field, we must make in the ground 
state energy of an unrestricted system, given in Eq. the substitution P^ — P^T" + P^T^. 

The last important detail in our approach refers to the orbital basis {4> v \v = l...k} used in our calculations. Instead 
of the spherical Bessel functions of Eq. 1)27(1 . the radial part of each orbital is decomposed in a sum involving five 
Gaussians confined to a sphere of radius Ro, while the angular part is maintained as defined by its symmetry. So, we 
change the basis in Eq. (|27|) by 

n— 1 5 

Ki,mi(*, ^ = N m (1 - x l II &i - x) £ V k e- D * R2 ° x2 Y hmi (6, 0), (30) 

i=l k=l 

where N n i is the orbital normalization, the polynomial in (1 — x) satisfies the boundary condition at x = 1 (r = Ro), 
and the polynomial in x makes the functions having / > be zero at the origin x = 0; also, the product in (an — x) 
makes the function be zero at the zeros an of the respective spherical Bessel function transposed to the interval 
< x < 1, and the last sum involves an expansion in five Gaussians. For N < 40, higher number of Gaussians in 
the expansion does not show any improvement on our results. The coefficients Vk as well as the exponents Dk are 
determined for each value of Ro, and are obtained by using the condition of maximizing the superposition between 
Eq. (|30fl and the respective spherical Bessel function. Once 14 and are determined and the basis is found, we run 
our RHF (UHF) code for given values of Ro and N, and find the parameters C V i (C"/' 5 ) that better describe Eq. © 
(Eq. QUI) and give the minimal energy in Eq. l(T3|) (Eq. 
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Finally, we would like to calculate two closely related quantities that give interesting information on the charging 
processes of a confined system. The first one is the QD chemical potential, fj, dot (N), which yields the energy difference 
between two successive ground states, and can be calculated as 

V dot (N)=E (N)-E (N-l). (31) 

The second one is the QD charging energy, E c h ar (N), which yields the energy cost for the addition of an extra electron 
to the system, 

E char (N) = I(N) - A(N) = E Q (N + 1) - 2E Q (N) + E (N - 1). (32) 

Here, I(N) = E (N - 1) - E (N) is the ionization potential, while A(N) = E (N) - E (N + 1) is the electronic 
affinity. It becomes clear, from these two equations, that E c har{N) = H dot {N + 1) — fi dot (N). 

IV. RESULTS 

Figure shows a comparison between the exact orbitals described by the spherical Bessel functions of Eq. 127|) 
and the expansions involving the Gaussians of Eq. H30JI used in our calculations for a GaAs QD having N — 2 and 
Rq = 100 A. The right panel shows the radial wavefunctions while the left one shows the radial probability densities. 
Although one can observe a difference in the s-wavefunctions, near x = 0, their probability densities are reasonable 
in that region. For any other occupation and radius, as also for InSb QDs, this same feature regarding s-orbitals is 
observed. Table [i] shows the optimized coefficients and exponents for the five Gaussians related to the six orbitals 
taken into account. The factor 10~ 4 in all exponents cancels the term Rq in Eq. (|30|l . Besides, all exponents related 
to any orbital having n — 2 are negative since there are regions where those wavefunctions possess negative values. 

In figure [21 we show results of a RHF-Roothaan calculation. In the left panel we present the ground state energy 
as function of the radius for a GaAs QD populated with 2, 8 and 18 electrons, while the right one shows the same, 
but for a InSb QD populated with 2, 34 and 40 electrons. These five distinct values, plus N = 20, are the successive 
magic numbers for the closed shell configurations at this symmetry. In both QDs, it is noticed that the kinetic energy 
is much higher and totally dominates the Coulomb one at smaller radii. Thus, the energies related to different N 
values are more distant from each other. At larger radii, or smaller electronic density, the electron-electron interaction 
becomes more important and the energy separation decreases. In the insets, which show a zoom at larger radii, we 
analyze the influence of the magnetic field on the system shell configuration. For both occupations N, the curves 
for increasing energy ordering refer to fields of 0, 2, 5, 8 and 10 T, even though we have only labeled the cases with 
N = 18 and N = 40. For N = 2 the presence of the magnetic field is imperceptible at those inset scales. As also 
expected, due to its large (/-factor, the Zeeman splitting is much higher for the InSb QD (observe the different energy 
scales in the insets). Most interesting is the fact that even for a closed shell configuration, where L = S = 0, the 
influence of the quadratic field (diamagnetic term) of Eq. I|26l which is the only non-zero contribution in the restricted 
case, becomes important mainly for higher occupation numbers and fields, as well as at larger radii (notice that there 
are almost no difference between and 2 T for any N). 

In order to study spectra for any occupation we show, in figure |3 the results of a UHF-Pople-Nesbet calculation for 
a GaAs QD having Rq — 100 A and without magnetic field. In the left upper panel it is made a comparison between 
the UHF results and the non-interacting electron case, where the energy shell structure is clear for N = 2, 8, 18, 20, 
34 and 40. Notice that the electron-electron interaction makes the energy of a QD, whose occupation corresponds 
to a shell less (more) than half-filled, be decreased (increased) with respect to the non-interacting case. When such 
occupation corresponds exactly to half-filled cases (N = 5, 13, 27 and 37), the interacting energy is approximately 
equal to the non-interacting one. 

In the left bottom panel of figure[3]we show both QD chemical potential (left scale, Eq. (^lfl ) and charging energy 
(right scale, Eq. (|32|l 'l. where the respective values of Eq are obtained from the unrestricted calculation presented in 
the left upper panel. Notice that n dot increases linearly as the occupation gets higher inside a given shell. When such 
shell is totally filled, there is an abrupt change in fi dot indicating that the following shell starts its occupation; observe 
that the higher the occupation, the larger is the changing. An anomalous behavior seems to occur for the 2s shell, 
whose n dot value is larger than the one for the If shell (that has higher energy). The charging energy is another way 
to verify not only the presence of shell structure for the spectrum, but also the validity of Hund's rule for the filling of 
such shells. In principle, E c h ar must present larger (smaller) peaks when the total (half) occupation of a given shell 
is achieved; the first fact is due to the higher difficulty to the addition of an extra electron to a QD when a filled shell 
state is reached. The second one refers to Hund's rule, which establishes that electrons must be added to the system 
with their spins being parallel, until all possible orbitals inside a given shell be occupied; this makes the total energy 
of the system smaller since this procedure maximizes the negative exchange contribution. However, some violations 
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FIG. 2: Restricted ground state energies for a GaAs (left panel) and InSb (right panel) QD populated with N = 2, 8 and 
18 (N = 2, 34 and 40) electrons. The Coulomb contribution is more important at larger radii. In the insets we analyze the 
influence of a magnetic field in the system, where one can see that even in a closed shell configuration having L = S = 0, the 
presence of the field is visible due to the diamagnetic term. 



of this rule can be verified in E c har- the smaller peak of N = 27 occurs here at N = 26, and the larger peak of N = 20 
has a negative value. 

The bottom and upper right panels of figure[3]show respectively the evolution of the total spin S and its projection 
Ms as a function of QD population, calculated from Eqs. (|23(l and (I24|) for the unrestricted energies. Notice that, 
without magnetic held, the Hund rule seems to be followed for all states with up to 40 electrons. The Ms expected 
value oscillates from in a hlled shell to its maximum in a half-filled shell, when it starts to decrease again on the 
way to the closing of the shell; the maxima are Ms — 1/2, 3/2, 5/2 and 7/2 for s, p, d and / shells, respectively. 
The S expected value yielded by the unrestricted formalism is also very reasonable; discrepancies are only observed 
at N = 24, where S > 2, and at N = 21, where S > 1/2. We believe that both discrepancies related to the 2s shell 
or to its surroundings - n dot larger than the one of 1/ shell, negative peak for N — 20 in E c h a r, and almost doubled 
S expected value for N = 21 - are caused by the non-reasonable Gaussian reproduction of this orbital, as visible in 
figure n These same qualitative results are observed for an InSb QD without magnetic field. 

By focusing on the lp shell we show in figurefl] for the same QD of the previous figure, how a finite magnetic field 
is able to violate Hund's rule in the system. Panels from left to right and from up to bottom show the successive 
ground state energies from N = 3 to N = 8 as this shell is filled, always considering that the Is shell remain fully 
occupied by two electrons, one spin-up and one spin-down; the distinct possible spin configurations for each TV are 
indicated by + (spin- up) and — (spin-down). In addition to the small Zeeman effect present in all occupations, there 
is a change of ground state spins at N = 4, 5 and 6 as the field is increased. Notice that at zero field the spin sequence 
is 1/2-1-3/2-1-1/2-0; in a field above 3 T it becomes 1/2-0-1/2-0-1/2-0, meaning that quartets 
and triplets are suppressed by the magnetic field, and the ground state of the system starts to oscillate only between 
singlets and doublets at high fields as N increases. When this lp shell is half- filled (N = 5), the ground state goes 
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FIG. 3: Unrestricted ground state energies for a Ro = 100 A GaAs QD without magnetic field. In the upper left panel we 
compare the unrestricted and the non-interacting energies, where the QD energy shell structure is visible. The bottom left 
panel shows QD chemical potential (left scale) and charging energy (right scale); the former displays abrupt change always that 
a new shell starts to be populated, while the latter presents larger (smaller) peaks when a shell is filled (half- filled), a direct 
consequence of the Hund rule. The bottom and upper right panels show, respectively, the ./V-evolution of the expected values 
of total spin S and its projection Ms- 



from a quartet to a doublet at Bq ~ 2 T; when it has one electron more (N = 6) or less (N = 4) than that, it goes 
from a triplet to a singlet at Bo ~ 3 T. 

At last, we have compared the results from both RHF and UHF self-consistent matrix formulations with the ones 
obtained from the LS-coupling scheme used in Ref. (3f|, where a GaAs QD having Rq = 90 A was considered, and 
the quadratic term in Bq was neglected since only small magnetic fields were considered. Also, only N — 2 and N = 3 
occupations were calculated, since the states were constructed analytically (not only a single Slater determinant), and 
the electron-electron interaction was included by using perturbation theory, that is justified at such radius. At zero 
field the energies for N = 2 are 16.5 meV (LS) and 16.1 meV (RHF), while for N = 3 they are 34.8 meV (LS) and 
33.9 meV (UHF). Therefore, both Roothaan and Pople-Nesbet formalisms indeed give smaller ground state energies 
than the LS perturbative scheme. We have also checked the validity of neglecting the diamagnetic term (quadratic in 
Bq) for fields smaller than 2 T. We may emphasize here that a disadvantage of the UHF approach is, in principle, that 
one is never sure to get trustable information about the expected values for L and S in a given QD state. On the other 
hand, the applicability of the LS scheme is cumbersome and becomes very complicated to be handled analytically as 
the QD occupation increases. 



V. CONCLUSIONS 



We have shown how the mean-field Roothaan and Pople-Nesbet formalisms applied to a spherical quantum dot 
confined system under applied magnetic field yield a fairly good description of its energy shell structure. For a 
maximum population of 40 electrons considered, the appropriated Gaussian basis set for each radius has been found. 
We have seen how a magnetic field influences the total energy of ground states even in closed shell configurations. We 
have also shown how both chemical potential and charging energy reproduce the closing and half-closing structures 
of the quantum dot energy shells. With the calculation of the total spin expected value for each occupation, in a 
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FIG. 4: Violation of Hund's rule induced by magnetic field in the Ro = 100 A GaAs QD of the previous figure. The panels 
show the successive occupation (indicated in the upper right corner of each panel) of the lp shell, assuming that the Is shell 
remains populated by one spin-up and one spin-down electron. The possible spin configurations for given N are indicated by 
+ (spin-up) and — (spin-down). For Bo = the spin sequence is 1/2 — 1 — 3/2 — 1 — 1/2 — 0, while at fields higher than 3 T it 
changes to 1/2 — — 1/2 — — 1/2 — 0. It is interesting to observe that a magnetic field is able to suppress the energy triplets 
and quartets of p-shells from the QD spectrum. 



given radius, we have observed that Hund's rule is satisfied at zero field. However, under finite magnetic field, we 
have shown that its applicability is violated and, at given values of the field, which depend on quantum dot radius 
and material parameters, there are transitions that change a given ground state symmetry. 
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